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Abstract 

An important question is if the gap in the high temperature cuprates has 
d x 2_ y 2 symmetry, what does that tell us about the underlying interaction re- 
sponsible for pairing. Here we explore this by determining how three different 
types of electron-phonon interactions affect the d x 2_ y 2 pairing found within 
an RPA treatment of the 2D Hubbard model. These results imply that in- 
teractions which become more positive as the momentum transfer increases 
favor d x 2_ y 2 pairing in a nearly half-filled band. 
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Typeset using REVTgX 



There have been a great deal of interest in the interpretation of recent experiments 
which address the question of the symmetry of the gap in the high temperature cuprate 



this would imply about the pairing mechanism. Here we discuss the relationship of the 
d x 2_ y 2 symmetry to the pairing mechanism by examining how three types of electron-phonon 
interactions affect d x 2_ y 2 pairing. 

As discussed by a number of authors a d x 2_ y 2 gap naturally arises from the exchange 
of antiferromagnetic spin fluctuations. However, the physical picture that emerges from these 
calculations is more general and shows that a d x 2_ y 2 gap will occur for a nearly half-filled 
band when there is an effective singlet interaction which is repulsive for on-site pairing and 
attractive for near-neighbor pairing. This spatial structure of the interaction means that its 
Fourier transform becomes more positive as the momentum transfer increases towards large 
values. This is easily understood from the BCS gap equation 



Near half-filling of the 2D system, the phase space is such that the important scattering 
process take electrons from (p f, — p 1) with p near a corner of the fermi surface, say 
near (tt, 0), to (p' |, — p' 1) with p' near (0, tt) or (0, — 7r). Since the interaction V(p — p') 
is positive, the relative phase of the states (p f , — P I) making up the bound Cooper pair 
changes sign as p goes from (ir, 0) to (0, tt) or (0, — tt) leading to a gap with d x 2_ y 2 symmetry. 
This is the case within an RPA approximation in which the interaction is mediated by the 
exchange of antiferromagnetic spin-fluctuations. It has also been found by Monte Carlo 
calculations |4]] that for the Hubbard model the pairing interaction is attractive in the d x 2_ y 2 
channel. In these cases the interaction is positive at all momentum transfers becoming larger 
in the region near (tt, tt) associated with the short range antiferromagnetic correlations. 

In order to explore the effect of electron-phonon interactions, we begin with a 2D Hub- 
bard model on a square lattice. 



superconductors []I|-£§. If this symmetry turns out to be d x 2_ y 2, it is natural to ask what 
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H = -t (4sCjs + c) s c is ) +UJ2 n^n il} (2) 
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Here t is a near neighbor hopping and U is the onsite Coulomb interaction. We take a 
simple phenomenological RPA form |J of the singlet pairing interaction associated with the 
exchange of spin fluctuations 

VsAp I -p) = \u\{ P '-p) (3) 

with x(q) — Xo(<z)/[l — Uxo{q)]- Here p = (p,iu n ), U is a renormalized Coulomb interaction 
and xo is the spin susceptibility 

/ \ 1 /( £ P+q) ~~ f( £ p) f A \ 

Xo(q, w) = — 2^ p (4) 

with e p = —2t(cosp x + cospy) — fi. Now, in addition to Vsf, we will examine three model 
electron-phonon interactions. The first is a simple on-site Holstein coupling of the form 

V x = 9Xirii (5) 

i 

with Xi the atomic displacement at site % and n% = + the onsite electron density. One 
could imagine this type of coupling arising from the interaction with an apical oxygen 0(4). 
In Eq. (5) the coupling is linear in the atomic displacement rather than quadratic. This is 
possible since 0(4) breaks the reflection symmetry with respect to a single Cu02 layer. The 
second electron-phonon interaction can be viewed as arising from the in-plane breathing 
motion [10| of an 0(2) oxygen 

v 2 = J29l x i( n i - n i+x) + Vi{ni - n i+y )). (6) 

i 

Here Xi describes the displacement of the 0(2) along the x-axis between the Cu sites at i 
and i + x and yi the y-axis displacement of an 0(2) along the y-axis between the Cu sites 
at i and i + y. The third interaction involves an axial z-motion of a buckled 0(2) atom. 

V 3 = J29l z U n i + n i+x) + 4i n i + n i+y)]- ( 7 ) 
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Here zf is for an 0(2) between the i and i+x sites and z\ is for an 0(2) between the i and i+y 
sites. Linear coupling of this type is possible for buckled Cu-O-Cu bonds. The electron- 
phonon interactions considered here are diagonal in the electron number representation. 
We note that it would also be interesting to study the effects of phonon modes where the 
electron-phonon coupling is not diagonal in the electron number representation. 

Assuming for the discussion that the lattice coordinate can be described as a local har- 
monic oscillator with frequency u (which just as g is of course different for the differ- 
ent modes), the effective electron-electron interaction mediated by the exchange of these 
phonons is 



V = TV n 7 (8) 

< + U 2 



where uj m is the Matsubara frequency 2mnT. Here for the local interaction Eq. (5) 
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^ 2 = ik (9) 



while for the breathing mode, Eq. (6), 



9^)\ 2 = ^-(sm 2 q x /2 + sm 2 q y /2) (10) 



2Mu 

and for the axial mode, Eq. (7), 

Kq)| 2 = ^£- Q (cos 2 q x /2 + C os 2 q y /2) (11) 

with M the O ion mass. 

In order to see how these interactions affect the pairing we examine the leading eigenvalue 
A and eigenfunction <p(p) of the Bethe-Salpeter equation, neglecting self-energy contribu- 
tions 0, 

\m = ~ £ [v SF ( P - pO + v ph ( P - p'WM-p'Mp'), (12) 
p' 

where G(p) is the single-particle Green's function given by 

G(p) = . (13) 

%Ul n Bp 



Here and in the following we will measure energies in units of t. The chemical potential has 
been chosen so that the site occupation (n^ + n^) = 0.875 and an effective Coulomb inter- 
action U = 2 has been taken. We will also take ujq = 0.25 and \g\ 2 /Mujq = 1 corresponding 
to an electron-phonon coupling strength \g\ 2 N(0)/Mu)Q ~ 0.4, where N(0) is the electron 
density of states at Hf- 

We find that the leading eigenvalue in the even frequency singlet channel has d x 2_ y 2 
symmetry and the temperature dependence of the eigenvalue \ x 2_ y 2{T) is shown in Figures 
1 and 2 for the electron-phonon interactions given by Eqs. (6) and (7), respectively. The solid 
line in each figure shows the eigenvalue in the absence of the phonon mediated interaction 
((7 = 0), while the dashed curve shows the effect of including the phonon mediated term. It 
is clear that the breathing mode interaction, Eqs. (6) and (10), suppress d x 2__ y 2 pairing 
while the axial 0(2) mode of Eqs. (7) and (11) enhances the d x 2_ y 2 pairing, raising T c . The 
local interaction, Eqs. (5) and (9), is orthogonal to the d x 2_ y 2 gap and hence does not affect 
the d x 2_ y 2 eigenvalue when self-energy effects are neglected. Including it in the self-energy 
will act to suppress T c due to the wave function renormalization |7]]. To understand the 
behavior shown in Figs. 1 and 2 we note that the strength of the coupling to the axial mode, 
Eq.(ll), decreases as q approaches (ir, ir). Because the phonon mediated interaction, Eq.(7), 
is negative, decreasing the magnitude of the coupling |g(q)| 2 acts to make the interaction 
more positive as the momentum transfer increases. As discussed in the introduction, this 
is the criteria for a d x 2_ y 2 gap to form when the system is near half-filling. Clearly it will 
be interesting to examine the isotope effect within models in which the strength of the 
electron-phonon coupling decreases at large momentum transfers. 

Thus we conclude that a d x 2_ y 2 gap implies that the pairing interaction becomes more 
positive for large momentum transfers. This is clearly the case for the spin-fluctuation 
interaction, Eq.(2), but as shown it can also occur for the attractive phonon mediated 
interaction if |<?(q)| 2 decreases at large momentum transfers. This form of coupling would 
also give rise to an electron-phonon coupling constant A which could be large compared to 
the effective coupling constant A tr entering transport processes since the transport coupling 
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constant At r weights large momenta transfers more heavily |HJ. For small momentum 
transfers q, a scattering of (p |, — p I) to (p + q T) — P — Ql) with p near a corner of the 
Fermi surface connects regions which have the same sign of the d x 2_ y 2 gap so that according 
to Eq. (1) an attractive electron-phonon interaction (negative V(q)) enhances A p . Another 
way to see that this latter case is similar to the spin-fluctuation interaction is to add U 
onto the phonon interaction giving U + Vph- For the d x 2_ y 2 channel, a constant has no 
effect but if it is larger than the magnitude of the phonon mediated interaction, Eq.(7), then 
as | C<3.) | 2 decreases, the total interaction U + Vph is positive and increases as q becomes 
large, just as Vsp. In order to obtain more quantitative information on the role of the 



electron-phonon interaction, it would be useful to have band structure calculations [jTT[] of 
the d x 2_ y 2 electron-phonon coupling constant for the i/-mode 

2E fe)fc ' g(k)g(k')^5(e k )5(e k ,) 

^g 2 {k)8{e k ) • 1 ) 

Here g(k) ~ (cosk x — cosk y ), e k is the band energy with k including the band index and 



M% k , is the electron-phonon matrix element for the z/'th phonon mode [12 
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FIG. 1. The d x 2_ y 2 eigenvalue of Eq. (2) versus the temperature T in units of the hopping t. 
The solid curve gives the eigenvalue for just the spin-fluctuation interaction Vsf and the dashed 
curve shows the effect when the electron-phonon interaction V2 associated with the breathing 
mode, Eqs. (6) and (10), is added to Vsf- The coupling constants are given in the text. 





